Effects of photon losses on phase estimation near the Heisenberg limit 
using coherent light and squeezed vacuum 
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Two path interferometry with coherent states and squeezed vacuum can achieve phase sensitivities 
close to the Heisenberg limit when the average photon number of the squeezed vacuum is close to 
the average photon number of the coherent light. Here, we investigate the phase sensitivity of such 
states in the presence of photon losses. It is shown that the Cramer-Rao bound of phase sensitivity 
can be achieved experimentally by using a weak local oscillator and photon counting in the output. 
The phase sensitivity is then given by the Fisher information F of the state. In the limit of high 
squeezing, the ratio (F — N) /N 2 of Fisher information above shot noise to the square of the average 
photon number N depends only on the average number of photons lost, ni oss , and the fraction of 
squeezed vacuum photons fi. For /i = 1/2, the effect of losses is given by (F — N)/N 2 = l/(l+2ni oss ). 
The possibility of increasing the robustness against losses by lowering the squeezing fraction /i is 
considered and an optimized result is derived. However, the improvements are rather small, with a 
maximal improvement by a factor of two at high losses. 

PACS numbers: 42.50.St, 03.65.Ta, 06.20.Dk, 42.50.Dv 



I. INTRODUCTION 

Quantum states of light can improve the sensitivity 
of phase measurements beyond the limits that apply to 
classical light sources. The phase sensitivity of coher- 
ent light (and hence of all classical light) is limited by 
the shot noise of independent photon detection events 
to the standard quantum limit of 5<fr 2 = 1/N. This 
limit can be overcome by using the multi-photon co- 
herences of non-classical light 0, H i, 0, [111, 0, H, !, 
[ToL ITU H3, EH. For two mode iV-photon systems, the 
highest possible phase sensitivity is achieved by max- 
imally path-entangled states, which are superposition 
states (|JV; 0) + \0;N))/y/2 where all photons are ei- 
ther in one path or in the other path of a two path 
interferometer HI QJ, GJ, [H GB [3, M, SH, HI- The 
phase sensitivity of these states defines the Heisenberg 
limit of 8(f) 2 = 1/N 2 . Since no TV-photon states achieve 
a higher phase sensitivity, this is the absolute limit of 
phase estimation for a fixed photon number N (Tl| . 

Unfortunately, it is rather difficult to generate maxi- 
mally path-entangled states u sing the available sources 
of non-classical light [23|, [HI HSl- It was therefore a 
significant discovery that the interference of a coherent 
state and a squeezed vacuum produces a high fraction of 
maximal path-entangled states when the average photon 
number from the squeezed vacuum is about equal to the 
average photon number of the coherent light [2(| [27| • In 
particular, Pezze and Smerzi showed that conventional 
two-path interferometry can achieve phase sensitivities 
close to the Heisenberg limit even in the presence of fluc- 
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tuating total photon number [27j ■ These results seem to 
put Heisenberg limited phase estimation within the reach 
of well-established quantum technologies. However, max- 
imally path-entangled states are very sensitive to photon 
losses, since the loss of a single photon can completely 
randomize the multi-photon coherence between the paths 
[H, US Hil . We can therefore expect that the Heisenberg- 
limited phase sensitivity achieved by coherent light and 
squeezed vacuum will rapidly decline as photon losses in- 
crease. 

In the following, we investigate the effect of photon 
losses on the phase sensitivity of the two-mode states 
generated by interference of coherent light and squeezed 
vacuum in detail. Assuming equal losses in both optical 
modes, we derive the mixed state after losses and find an 
optimal phase estimator based on the general analysis 
for quantum metrology using mixed states [3l[. We find 
that the Cramer-Rao bound giving the maximal phase 
sensitivity of the state can be achieved by a simple ex- 
perimental setup using a weak local oscillator field and 
photon detection. It is therefore possible to obtain a 
phase sensitivity equal to the Fisher information F in 
an experimentally feasible setup using only linear optics 
and photon detection. In the limit of high squeezing, 
the ratio (F — N) /N 2 of Fisher information above shot 
noise to the square of the average photon number N de- 
pends only on the fraction of squeezed vacuum photons 
fi and the average number of photons lost ni oss . Thus the 
effect of photon losses on phase sensitivities close to the 
Heisenberg limit has the same dependence on the average 
number of photons lost, regardless of the total average 
photon number N. In particular, photon losses reduce 
the phase sensitivity for Heisenberg-limited estimation 
at fi — 1/2 by a factor of 1/(1 + 2ni oss ). Finally, we 
investigate the possibility of improving the Fisher infor- 
mation by optimizing the squeezing fraction /i for a given 
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FIG. 1: Illustration of Heisenberg limited estimation of a 
small phase shift cj> with coherent light and squeezed vacuum 
in the presence of linear losses. The probability that any given 
photon is lost between generation and detection is given by 
the loss fraction a. 



number of photons lost. However, the result shows only 
small improvements, approaching a maximal increase of 
Fisher information by a factor of two in the limit of high 
photon losses. 



II. EFFECTS OF LOSSES ON THE TWO MODE 
SQUEEZED-COHERENT STATE 

Fig. [1] shows a possible experimental setup realizing 
Heisenberg limited phase estimation with coherent light 
and squeezed vacuum in the input ports of a two-path 
interferometer. Initially, mode d\ is in a coherent state 
|a) and mode 02 is in a squeezed vacuum state \S(r)) — 
S r (r)|vac), where S(r) = exp[l/2(ra2a2 — r*a 2 a 2 )} is 
the squeezing operator. Interference at the initial beam 
splitter then results in multi-photon coherences between 
the two paths inside the interferometer, as discussed in 
[H, H3- However, photon losses occurring at any point 
between the generation and the detection of the light 
fields will reduce these multi-photon coherences. In the 
following, we assume linear losses with equal loss rates in 
the two modes. It is then possible to represent the losses 
by the loss fraction <r, defined as the probability that any 
given photon is lost between generation and detection. 

The effects of linear losses on the two input modes 
correspond to interference with a vacuum state, followed 
by a trace over the modes representing the losses. Since 
photon losses from orthogonal modes are statistically in- 
dependent, it is possible to consider the effect of photon 
losses on the two input modes separately. For the co- 
herent state, the losses simply reduce the amplitude a 
by a factor of \/T~ 

Ored. 



a, so that the output amplitude is 
VT — era and the density matrix pi of mode Si 
after losses is 

Pi = |a r cd.)(arcd.|- (1) 
In the case of the squeezed vacuum, losses change the 



variances of the quadrature components &i and j/2 of the 
field mode hi = £2 + ijj2- The output is a Gaussian state 
with quadrature variances of 



<J + (1 
4Ay 2 2 = a + 



a)e~ 2r 
a)e 2r . 



(2) 



In general, a Gaussian mixed state defined by the vari- 
ances Ax 2 , and Ay 2 can be described by a squeezed ther- 
mal state, 



f>2 = S'(?'rcd.)Pth(A)S' t (r rcd .), 

where the thermal state is given by 

00 

/5 th (A) = (l-A)£A»(n| 



(3) 



(4) 



and r ro d. < r is the reduced squeezing parameter ob- 
tained from the ratio of the variances after losses. In 
terms of the parameters A and r ro d., the variances of pi 



are 



<±Axl 



4Ay 2 2 = 



1 + A 

T-a ( 

1 + A 



(5) 



Thus it is possible to determine the values of A and r re d. 
corresponding to an initial squeezing parameter r and a 
loss probability a from eqs. ([2]) and ([5]). 

The complete two mode state after losses is given by 
the product of the states in mode a\ and in mode a 2 , 



p = P\® P2- 



(6) 



The phase sensitivity achieved by this state can be ana- 
lyzed using the general formalism for mixed states [3l[ . It 
is then possible to determine both the quantum Cramer- 
Rao bound of phase estimation and the measurement 
procedure that achieves this bound in the presence of 
photon losses. 



III. PHASE ESTIMATION WITH MIXED 
STATES 

Quantum phase estimation is performed by measuring 
a phase estimator A. If the average value of A is chosen to 
be zero at <fi = 0, the phase derivative of the average gives 
the ratio of the average of A and the small phase shift 
that quantum estimation seeks to detect. Thus the av- 
erage phase estimate (0 cs t.) = {A) /{d{A) /dej)) converges 
on the correct value <fi as the number of measurements 
increases. However, each individual measurement has a 
statistical error of 8(f> 2 that determines how quickly the 
average of the measurement results converges on the cor- 
rect value of d>. In terms of this measurement error, the 
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phase sensitivity obtained with a specific estimator A is 
given by 



Tr 



(7) 



An optimal phase estimator A maximizes this phase sen- 
sitivity and achieves the quantum Cramer-Rao bound of 
the state p. As was shown in 31| , the optimal estimator 
is given by the symmetric logarithmic derivative G of the 
density matrix p, as defined by the operator relation 



d 



Gp 



(8) 



presence of losses. The density matrix was derived in sec. 
HI! and the generator h for a two path interferometer is 
given by half the photon number difference between the 
two paths. In the present context, the two-mode density 
matrix is a product of the two input mode density matri- 
ces. It is therefore convenient to express the generator h 
in terms of the input modes Si and 0,2, which are equal 
superpositions of the modes describing the two paths in- 
side the interferometer. The photon number difference 
between the two paths is then equal to an interference 
term of the input modes a\ and 0,2. Specifically, it can 
be written as 



a 2 d\ 



(12) 



Note that this relation does not uniquely define G if p has 
eigenvalues of zero. In that case, any operator G fulfilling 
eq. (J8|) is an optimal estimator. The maximal phase 
sensitivity achieved by an optimal estimator G is equal 
to the Fisher information F = 1 /6^r pt of the quantum 
state p. The Fisher information can be evaluated from 
eq. (J7|) by using A — G and eq. ([8|). The result is equal 
to the variance of the optimal phase estimator G, 



F = Tr 



{#>} 



(9) 



Eqs. © and © summarize the results for quantum 
metrology with mixed states obtained in [3l[ without the 
explicit expansion into eigenstates of the density matrix 
used in the original derivation. In general, these results 
apply to any parameter (f> that changes the quantum state 
p. In the specific case of a phase shift, <j) is the parameter 
of a unitary transformation exp[— icfih] generated by an 
operator h. The phase derivative of the density matrix is 
therefore given by the commutation relation of p and h, 



d_ 
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p = -i 



(hp — ph 



(10) 



To find an optimal estimator G for a given generator 
h and a given quantum state p, we have to solve eq. 
© using the phase derivative given by eq. (JTDJ) . This 
relation can be summarized by 
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For unitary transforms, an optimal estimator G is there- 
fore obtained when the anti-commutation of G and p has 
the same form as the commutation of h and p. 



IV. DERIVATION OF AN OPTIMAL PHASE 
ESTIMATOR 

We can now derive an optimal estimator for phase es- 
timation with coherent light and squeezed vacuum in the 



With this generator, eq. (jTTJ) provides a relation between 
the optimal estimator G and the two mode state p in 
terms of the creation and annihilation operators of the 
input states, 



Gp + pG 



(a\a 2 



a}a 2 



alai 



p. (13) 



Since the density matrix can be written as a product of 
states in mode a\ and mode 02, the effects of the opera- 
tors eti and S2 on the states p\ and p~2 can be determined 
separately. It is then possible to derive a particularly 
simple form of G by only considering the relations be- 
tween single mode Gaussian states and the creation and 
annihilation operators of their respective modes. 

First, we consider the effects of the annihilation and 
creation operators of mode d\ on the coherent state den- 
sity matrix p\. The coherent state | a re d.) is a right 
eigenstate of the annihilation operator a\. It is therefore 
possible to replace the operator d\ operating from the 
left on p and the operator a\ operating from the right 
with the complex number a 10 d.- The application of a\ 
to the coherent state | a re d.) changes that state into a 
superposition of the original state with an amplitude of 
a re d. and an orthogonal state that can be represented by 
a displaced one photon state. The creation operator can 
thus be written as a sum of the coherent amplitude a re d. 
and an operator that changes the initial state into an 
orthogonal state, 



it 



ared 



a Ie d. 



(14) 



It is possible to separate the relation for the estimator G 
given by eq. (fT3| into two parts, one that leaves the co- 
herent state unchanged, and one that describes the tran- 
sition matrix elements between the coherent state and the 
displaced one photon state. The separation is achieved 
by writing the optimal estimator as G = g\ + 52 , where 
gi is the component of the estimator associated with the 
transition matrix elements in mode &i and §2 is the com- 
ponent of the estimator that commutes with the coherent 
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state pi. The two relations denning g± and c/ 2 then read 

hp + pgi = 

- ((a[ - a Ted .)a 2 p - pa\{a x - a rcd .)) (15) 



92 P + P92 



od. ((a 2 - a\)p - p(a 2 - 4)) 



(16) 



To find <?i, we make use of the fact that (di — a Ic d.)p = 
0. It is therefore possible to add or subtract multiples of 
this operator and its self-adjoint operator to the right side 
of eq. (fT5"|) without changing the relation. The solution 
for the self-adjoint operator g\ obtained in this manner 
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(17) 



To find §2 , we make use of the fact that eq. ([T6|) only 
includes operators acting on the state in mode d 2 . There- 
fore, we only need to consider the density matrix p 2 of 
the squeezed thermal state. Since (d 2 — d|) = 2iy 2 corre- 
sponds to the anti-squeezed quadrature component, the 
right side of eq. (flE|) corresponds to the commutation 
relation between the quadrature y 2 and the density ma- 
trix p2- For Gaussian states, this kind of commutation 
relation can be converted into an anti-commutation rela- 
tion for a different quadrature component. In the case of 
the squeezed thermal states given by p 2 , the conversion 
is given by 



i [hfc - P2V2) 



1 - A 
1 + A 



e 2r " d ■ {X2P2 + P2X2). (18) 



Comparison with eq. (| 16(1 indicates that g 2 is a multiple 
of the quadrature component i 2 . Specifically, eq. (lit 
can be solved by 
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92 = 2Axl X2 > 



(19) 



where the coefficients A and r re d. have been expressed in 
terms of the quadrature variance Ax| using eq. ([5]). 

The optimal estimator G is given by the sum of gi 
and §2 ■ Since <?i is a quadratic function of the creation 
and annihilation operators, and §2 is a linear function of 
the operators of mode d 2 , it is possible to express G as 
a quadratic function of the field operators. Specifically, 
the result can be written as an interference term of a 2 
and a field b, 



G = tf 



0.2 



where the field b is given by 
/ 1 



a re d. 



1 



ai. 



(20) 



(21) 



Experimentally, this estimator can be realized by sub- 
tracting a coherent amplitude of a re d.(l/4Ax2 + 1) from 
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FIG. 2: Schematic setup for the observation of optimal phase 
sensitivity obtained by coherent light and squeezed vacuum in 
the presence of losses. The optimal phase estimation is real- 
ized by subtracting a coherent amplitude of a rc d. (l/4Aa^2 + 1) 
from the output mode fii using a local oscillator and a beam 
splitter with low transmittivity T (BS3). 



the output field in di using interference with a local os- 
cillator. A possible setup is shown in figHJ At zero phase 
shift, the output modes correspond to the input modes. 
A small displacement of the coherent field mode di is 
realized by interference with a local oscillator field at a 
beam splitter of very low transmittivity T (BS3). Finally, 
the displaced mode b and the output mode a 2 interfere 
at a fourth beam splitter (BS4). The phase estimator is 
then equal to the photon number difference in the output. 
Interestingly, the estimator is a linear function of the de- 
tected photon numbers. This is quite different from the 
optimal phase estimation for pure states considered in the 
initial work on Heisenberg limited phase estimation with 
coherent and squeezed light, where higher order moments 
of the detected output photon number distribution were 
essential [l?} . The present setup therefore represents a 
major simplification of the phase estimation procedure 
for phase sensitivities close to the Heisenberg limit, even 
in the pure state case where the phase sensitivity is equal 
to that obtained from direct photon counting in the out- 
put. 



V. DEPENDENCE OF PHASE SENSITIVITY 
ON PHOTON LOSSES 

As mentioned in sec. IIII1 the Fisher information of 
the quantum state p is equal to the expectation value of 
the squared estimator G. Using the result of eq. (j!ifl)) . 
the Fisher information of the squeezed-coherent state p 
is found to be 



Tr{G 2 p} 



Tod. 



4 Ax 2 2 



ri2, 



(22) 



where n 2 is the average number of photons in the 
squeezed mode d 2 after losses. Here, the effects of losses 
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are expressed indirectly through the values of U2, cv rc( j. 
and Ax2- The specific effects of a loss probability of 
a on the input state are given by a ro( j. = \/l — era, 
ri2 = (1 — cr) sinh 2 r and eq. ©. The Fisher information 
can then be expressed in terms of the input amplitude 
a, the input squeezing r, and the loss probability a. The 
result reads 



F = (1 - a) 



(1 



-2r 



sinh 



(23) 



In this representation of the Fisher information, the most 
significant effect of the losses is the limitation of squeez- 
ing effects represented by e~ 2r . However, it is difficult 
to see how this limitation relates to the maximal phase 
sensitivities achieved at equal intensities of coherent light 
and squeezed vacuum. It is therefore convenient to ex- 
press the result in terms of average photon numbers in- 
stead. 

The total average photon number after losses is given 
by N = (1 — a)(a 2 + sinh 2 r). To evaluate the distri- 
bution of photons between the coherent light and the 
squeezed vacuum, we introduce the squeezing fraction 
fi = (1 — a) sinh 2 r/N, defined as the fraction of photons 
in the squeezed mode &2- Finally, the effects of losses 
can be given in terms of the average number of photons 
lost, ni oss = Na/(1 — a). Since N is the average photon 
number after losses, the total photon number of the in- 
put state is given by the sum of N and ni oss , as shown in 
fig. [3] The Fisher information is then given by 



F = N 



4(1 -A^ 



1 



-1r 



4fi n\ Q 



N. 



(24) 



Note that this phase sensitivity can be greater than N 2 , 
since the actual Heisenberg limit for fluctuating photon 
numbers is given by the average of the squared photon 
number, not the square of the average photon number 

Since we are mainly interested in Heisenberg limited 
phase sensitivities with large photon numbers, it is rea- 
sonable to assume that the squeezing levels will be high 
enough to satisfy e" 2r <C 1. We can then neglect the r- 
dependent term in eq. (|24p to obtain a particularly sim- 
ple relation between phase sensitivity and photon losses. 
Specifically, the Fisher information above the standard 
quantum limit, F — N, is given by a fraction of ./V 2 deter- 
mined only by the squeezing fraction /i and the average 
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FIG. 3: Illustration of the definition of TV and n\ 
average photon number after losses. 
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FIG. 4: Effects of photon losses ni oss on the enhancement 
of sensitivity (F — N)/N 2 for phase estimation at squeezing 
fractions of /i = 1/2 and n = 1/8, 1/32, 1/128. Since photon 
losses reduce the enhancement of sensitivity by a factor of 
1 + 4u nioss, states with lower squeezing fractions are more 
robust against photon losses than the states with maximal 
multi-photon coherence at /i = 1/2. 



number of photons lost. Since this fraction does not de- 
pend on N, it provides a photon number independent 
expression of the effects of losses on the phase sensitiv- 
ity of squeezed-coherent states. In the following, we will 
refer to this expression as the enhancement of sensitivity, 



F-N _ 4/i(l - n) 
N 2 ~ 1 + 4/i m ost 



(25) 



As reported in [26|, |27[, equal intensities of squeezed 
vacuum and coherent light (ji — 1/2) result in maxi- 
mal multi-photon coherences, including a significant frac- 
tion of maximally path-entangled states. In the absence 
of losses, the enhancement of sensitivity for these pure 
states is (F — N)/N 2 = 1, the maximal value that can 
be achieved in the limit of high squeezing. However, eq. 
<P5|) also shows that photon losses rapidly reduce this 
enhancement of phase sensitivity. Specifically, 



F-N 



N 2 



1 



1 + 2nio 



(26) 



This dependence of phase sensitivity on the average num- 
ber of photons lost reflects the fact that the loss of a single 
photon completely randomizes the iV-photon coherence 
of a maximally path entangled state, irrespective of the 
total photon number N. Thus, the average loss of just 
half a photon already reduces the enhancement of sensi- 
tivity to half of its original value. 

Eq. (|25p indicates that the effect of photon losses on 
the enhancement of sensitivity (F — N)/N 2 decreases 
when the squeezing fraction fi is lowered. Specifically, 
photon losses reduce the enhancement of sensitivity by a 
factor of 1 + 4/i ni oss , defined by the product of squeez- 
ing fraction and photon losses. Therefore, states with 
lower squeezing fraction [i are more robust against pho- 
ton losses. Fig. @] shows a comparison of the loss- 
dependent enhancements of sensitivity [F — N) /N 2 for 
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FIG. 5: Comparison of enhancement of sensitivity (F — N)/N 
obtained at the optimized squeezing fraction /i = /i op t. with 
the enhancement of sensitivity obtained at /i = 1/2. 



FIG. 6: Ratio of the sensitivity enhancement obtained at op- 
timal squeezing fraction and the sensitivity enhancement ob- 
tained at /i = 1/2. The ratio is close to one at low losses and 
approaches a maximal value of two at high losses. 



different squeezing fractions /i < 1/2. At low losses, 
the enhancement of sensitivity is maximal for fi = 1/2 
and decreases with decreasing /i. As losses increase, the 
enhancement of sensitivity for /i = 1/2 drops to values 
below the corresponding enhancements at lower /z, indi- 
cating that states with lower squeezing fraction can have 
higher Fisher information in the presence of losses. If 
the average number of photons lost is fixed, the highest 
enhancement of sensitivity given by eq. (|25|) is found at 
a squeezing fraction of 
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fJ-opt. 



4n 



loss 



(VI + 4n loss - 1) 



(27) 



The enhancement of sensitivity at this optimal squeezing 
fraction (J, op t. is given by 



N 



N 2 



\/l + 4ni oss - 1 

2^1oss 



(28) 



Fig. [5] shows a comparison of the enhancement of sen- 
sitivity at fi op t. with the enhancement of sensitivity at 
fx = 1/2. Although the reduction of squeezing fraction /i 
results in higher enhancements of phase sensitivity, the 
relative improvements seem to be rather small. Fig. [5] 
shows the improvement factor given by the ratio of the 
enhancement of sensitivity at /i op t. and the enhancement 
of sensitivity at /i = 1/2. The improvement factor is neg- 
ligibly small at low losses, with a value of only 1.072 at 
average losses of half a photon. Thus, the optimization of 
the squeezing fraction can do little to compensate the re- 
duction of the enhancement of phase sensitivity to half its 
value at ni oss = 1/2. As losses increase, the improvement 
achieved by an optimization of the squeezing fraction be- 
comes more significant. However, the improvement is 
limited by its asymptotic value of 



lim„ 



(F-N)\ 



Ai=Mo P t. 



(F-A% =1/2 



(29) 



so that the optimization of the squeezing fraction can 
at most double the enhancement of phase sensitivity 



achieved at a squeezing fraction of fj, = 1/2. The phase 
sensitivity achieved at [i = 1/2 therefore remains close to 
the maximal phase sensitivity that can be achieved with 
any squeezed-coherent state, even in the presence of very 
high photon losses. 



VI. CONCLUSIONS 

We have shown that photon losses reduce the phase 
sensitivity of the multi-photon coherences obtained from 
interferences of equal intensities of squeezed vacuum and 
coherent light by a factor of 1 + 2ni oss , where ni oss is the 
average number of photons lost. This result corresponds 
to the expectation that a single photon loss randomizes 
the coherence of maximally path entangled states, re- 
gardless of the total photon number N. A small improve- 
ment of the robustness against losses can be achieved by 
reducing the fraction of squeezed vacuum in the total 
photon number. However, the improvements are rather 
small and indicate that the robustness against losses of 
A-photon states depends mainly on their phase sensitiv- 
ity, regardless of the type of state used. 

We have also shown that the Cramer-Rao bound of 
squeezed-coherent states in the presence of losses can be 
achieved in experimentally feasible measurements using 
a weak local oscillator field, linear optics, and photon 
counting. Interestingly, the use of the local oscillator 
simplifies the phase estimation procedure to an estima- 
tor linear in the detected photon numbers. It may there- 
fore present an interesting alternative to direct photon 
counting in the output, even if the improvement in phase 
sensitivity is negligibly small. 

In general, our results confirm that phase estima- 
tion near the Heisenberg limit can be performed with 
squeezed-coherent light, but only if the average number 
of photons lost can be kept low. Therefore, high efficien- 
cies of photon transmission and detection will be essential 
for quantum metrology close to the Heisenberg limit. 
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